Microscopic self-energy of ^°Ca from the charge-dependent Bonn potential 
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The effects of short-range correlations on the nucleon self-energy in ''"Ca are investigated using the 
charge-dependent Bonn (CDBonn) interaction. Comparisons are made with recent results for the 
self-energy of '"'Ca derived from the dispersive optical-model (DOM). Particular emphasis is placed 
on the non-locality of the imaginary part of the microscopic self-energy which suggests that future 
DOM analyses should include this feature. In particular, data below the Fermi energy appear 
sensitive to the implied orbital angular momentum dependence of the self-energy. Quasiparticle 
properties obtained for the CDBonn interaction are substantially more mean-field-like than the 
corresponding DOM results with spectroscopic factors larger by about 0.2 e.g. Reaction cross 
sections obtained from the microscopic self-energy for scattering energies up to 100 MeV indicate 
that an adequate description of volume absorption is obtained while a considerable fraction of surface 
absorption is missing. The analysis of the non-locality of the imaginary part of the microscopic self- 
energy suggests that a simple gaussian provides an adequate description, albeit with rather large 
values for /3, the non-locality parameter. 



I. INTRODUCTION 

Calculations of single-particle (sp) observables in nu- 
clei starting from realistic nucleon-nucleon (NN) inter- 
actions have demonstrated that several energy scales in- 
fluence these quantities [Ij. Both short-range as well as 
tensor physics associated with energy scales far above and 
below the Fermi energy, in addition to long-range physics 
associated with the nuclear surface, and most influential 
near the Fermi energy, determine single-nucleon motion 
in the nucleus. This interpretation is well established 
for nuclei near stability where (e, e'p) experiments have 
yielded unambiguous results approaching absolute spec- 
troscopic factors [H. The observed reduction of the sp 
strength near the Femi energy of about 35% can only be 
quantitatively explained by a combination of the effects 
of long-range and short-range correlations (SRC) [H, 0]. 
While it appears that the bulk of the depletion is associ- 
ated with the coupling of sp motion to low-lying collec- 
tive states and giant resonances that are mostly associ- 
ated with surface phenomena, there is an approximately 
modest 10% depletion due to the short-range and tensor 
force of the NN interaction calculated from the micro- 
scopic self-energy obtained for ^^O These calcula- 
tions also demonstrate an accompanying modest appear- 
ance of high-momentum components in the ground state. 
The presence of about 10% high-momentum components 
does have a signiflcant influence on the binding energy 
of the ground state, in some cases accounting for about 
65% of the binding energy [5|. 

While the details of these high-momentum components 
depend sensitively on the chosen NN interaction, exper- 
imental evidence confirms their presence in about the 
amount predicted although the details of their en- 
ergy dependence require further study 0]. Calculations 
for "'^^O employing other methods that include the ef- 



fect of SRC confirm that they remove strength from va- 
lence orbits by about 10% [1, US]- The presence of high- 
momentum components in the ground state for nuclei 
heavier than ^^0 is usually obtained by performing a 
local-density approximation of the nuclear-matter spec- 
tral function obtained by the correlated-basis-function 
method at several densities [ill, ■ It appears that only 
the Green's function calculations of Refs. can gen- 

erate spectral functions directly for such nuclei. 

Green's function calculations generate self-energies 
that can also be employed at positive energy where 
they can be transformed into the so-called reducible 
self-energy that is equivalent with the elastic-nucleon- 
scattering T- matrix |13l | . Elastic- nucleon-scattcring data 
are usually analyzed with so-called optical potentials that 
provide quite accurate representations of the data [l^ 
llg . Optical potentials are normally formulated in terms 
of local potentials that have no immediate interpretation 
as nucleon self-energies. A close connection to the nu- 
cleon self-energy has been developed in the dispersive- 
optical-model (DOM) proposed by Mahaux and Sar- 
tor [itI . [Tsj . By employing the dispersion relation be- 
tween the real and imaginary part of the self-energy in 
a subtracted form, only the imaginary part of the self- 
energy and the real part at the subtraction point (typi- 
cally the Fermi energy) are required to provide a fit to the 
available data. These include not only elastic-scattering 
data but also bound-state information, thereby provid- 
ing a natural link between nuclear reactions and nuclear 
structure when appropriate assumptions are made about 
the relation of the imaginary part above and below the 
Fermi energy [l^ ■ Standard assumptions about the func- 
tional form of these potentials include surface and vol- 
ume contributions with conventional form factors. The 
parameters of these potentials are then constrained by 
a fit to the available data. Recent implementations of 
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this approach attempt a simultaneous fit to nuclei with 
different nucleon asymmetry thereby providing access to 
extrapolations to more exotic nuclei |19l - l21| . An exten- 
sion of the DOM was recently introduced which enhances 
the domain of its applicability for data below the Fermi 
energy [22|. Conventional DOM analyses transform the 
non-local real part of the self-energy at the Fermi energy 
into a local but energy -dep endent potential according to 
a standard procedure [2j|. Such an energy-dependence 
distorts the normalization of solutions of the Dyson equa- 
tion which is of particular relevance for the calculation of 
the spectral function below the Fermi energy. By return- 
ing to a non-local form and renormalizing the imaginary 
part in a well-defined manner it is possible to ob- 
tain properly normalized results for spectral functions, 
the one-body density matrix, etc. [2^ . 

Since all implementations of the DOM have assumed 
a local representation of the imaginary part of the self- 
energy, it is of particular relevance to investigate the non- 
local content of microscopically calculated self-energies 
that are based on a realistic interaction. A recent investi- 
gation of this typ e focused on the coupling of the nucleon 
at low energy [2J] . Such long-range correlations are well- 
described by the Faddeev-random-phase approximation 
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(FRPA) I, Hi-il]. The analysis of Ref. 
that at low energy, mostly representing the coupling to 
surface excitations, there is ample evidence of a substan- 
tial non-local contribution to the nucleon self-energy. In 
the present work we extend this investigation to study 
the predominantly volume coupling associated with the 
contribution of SRC as determined by the ab initio strat- 
egy of Refs. 

We employ the CDBonn interaction [sO, Ej in the 
present work which is a relatively soft interaction js^ . 
At the same time, we establish how such an ab initio 
optical potential fares in the description of elastic scat- 
tering. Since most recent DOM results include '^°Ca, we 
have also performed the present calculations for this nu- 
cleus. Wc note that fully microscopic calculations of the 
optical potential for finite nuclei are sorely lacking and 
most applications involve a detour through nuclear mat- 
ter accompanied by a suitable local-density approxima- 
tion jsll. The work reported in Ref. [l^ demonstrates 
that the DOM is capable of incorporating SRC in its 
description of the nucleon self-energy. Its current imple- 
mentation however, does not yet properly account for the 
effect that SRC have on the distribution of sp strength, in 
particular that high-momenta appear increasingly likely 
at more negative energies We therefore also con- 

sider the high-momentum components generated by the 
CDBonn interaction in "'"Ca to clarify this issue. 

In Sec. In] we introduce the essential steps involved in 
the microscopic calculation of the self-energy for ^°Ca. 
While the initial step involves a nuclear-matter calcula- 
tion, the subsequent steps correctly incorporate the prop- 
agation of nucleons in the finite system under considera- 
tion. The relevant quantities that can be derived by solv- 
ing the Dyson equation above and below the Fermi energy 



are also summarized there. Results for spectral functions, 
momentum distributions, natural orbits, charge density, 
energy sum rule, and elastic scattering are presented in 
Sec. mil Whenever relevant, these results are compared 
with corresponding DOM quantities. We also attempt 
to identify properties of the non-local imaginary part of 
the microscopic sclf-encrgy for future implementation of 
DOM analyses. Conclusions are presented in Sec. |IVl 



II. THEORY 

A. Effective interaction 

The microscopic calculation of the nucleon self-energy 
proceeds in two steps, as employed in Refs. We 
start with an outline of the general idea and initially 
employ a schematic formulation. A diagrammatic treat- 
ment of SRC always involves the summation of ladder 
diagrams. When only particle-particle (pp) intermediate 
states are included, the resulting effective interaction is 
the so-called C/-matrix. The corresponding calculation 
for a finite nucleus (FN) can be represented in operator 
form by 



GpNiE) + VGfj^{E)gFN{E), 



(1) 



where the noninteracting propagator G^jy (^) represents 
two particles above the Fermi sea of the finite nucleus 
taking into account the Pauli principle. The simplest 
implementation of G^'jy involves plane-wave intermedi- 
ate states (possibly orthogonalized to the bound states). 
Even such a simple assumption leads to a prohibitive cal- 
culation to solve Eq. ([T|) and subsequently generate the 
relevant real and imaginary part of the self-energy over 
a wide range of energies above and below the Fermi en- 
ergy. We are not aware of any attempt at such a direct 
solution at this time, except for the use of the ^/-matrix 
as an effective interaction at negative energy. Instead, 
we employ a strategy developed in Refs. 0, HH that 
first calculates a ^-matrix in nuclear matter at a fixed 
density and fixed energy according to 

Qnm{Emm) = V + VG^^j^j{Enm)Gnm{Enm)- (2) 

The energy E^m is chosen below twice the Fermi energy 
of nuclear matter for a kinetic energy sp spectrum and 
the resulting Qmm is therefore real. Formally solving 
Eq. ([l} in terms of Q nm can be accomplished by 



Gfn{E) — Gnm 

+ Qnm {G%{E) - G?,V} Qfn{E), 



(3) 



where the explicit reference to E^m is dropped. The 
main assumption to make the self-energy calculation 
manageable is to drop all terms higher than second order 
in f/ATA/, leading to 



Qfn{E) — Gnm ~ QnmG^j^j^jQnm 

-\- yNM Gjj,{E)gj,M, 



(4) 
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FIG. 1. Graphical representation of the BHF (a), the two- 
particle-one-hole contribution (b) and one-particle-two-hole 
term (c) to the self-energy of the nucleon. The 0jvA/-niatrix 
is indicated by the wiggly line. 



where the first two terms are energy- independent. Since 
a nuclcar-mattcr calculation already incorporates all the 
important effects associated with SRC, it is reasonable to 
assume that the lowest-order iteration of the difference 
propagator in Eq. ([4]) represents an accurate approxima- 
tion to the full result. This assertion does require further 
confirmation in future studies. 

The Bethe-Goldstone equation for Qnm reads in the 
appropriate basis 



{ki\ e^lf^^^ \k'e') = {k£\ 1^,^,^ 

+ \Y, I dk"{k"f (fc^ll/^-'^^^^l/fc'r') 

2 ,„ J 



(5) 
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The variables fc, A;', and k" denote the relative wave vec- 
tors between the two nucleons, I, t' , and £" the orbital 
angular momenta for the relative motion, K and L the 
corresponding quantum numbers for the center-of-mass 
motion, S and T the total spin and isospin, and Js is ob- 
tained by coupling the orbital angular momentum of the 
relative motion to the spin S. Wc note that Q{K, k) is 
the angle-averaged form of the product of step-functions 
that allows a partial wave decomposition. The implied 
sole dependence of the G-matrix on the magnitude of 
K , the center-of-mass wave vector, also ensures that the 
solution of Eq. ([5]) does not depend on L. The L-label 
is kept however, to facilitate the recoupling to individ- 
ual orbital angular momentum states, as discussed be- 
low. Equation ([5]) generates an appropriate solution of 
two-body short-range dynamics but the resulting matrix 
elements require further manipulation before becoming 
useful for the finite nucleus. 

The self-energy contribution of the lowest-order term 
Qnm in Eq. ([4]) is shown in Fig. [Ija) and is similar 
to a Brucckncr-Hartrce-Fock (BHF) self-energy. While 
strictly speaking the genuine BHF approach involves self- 
consistent sp wave functions, as in the HF approximation, 
the main features associated with using the t/jvj\/-matrix 
of Eq. ([S]) are approximately the same when employing a 
summation over the occupied harmonic oscillator states 
of ''^Ca. Hence we will use the BHF abbreviation. The 
correction term involving the second-order in Q nm calcu- 
lated in nuclear matter is also static and can be obtained 



from the second term in Eq. ([5]) by replacing the bare 
interaction by Qnm- The corresponding self-energy is 
also real and generated by summing over the occupied 
oscillator states in the same way as for the BHF term. 

The second-order term containing the correct energy- 
dependence for t/i?7v in Eq. (jl]) can now be used to con- 
struct the self-energy contribution depicted in Fig. [TJb), 
representing the coupling to two-particle-one-holc (2plh) 
states. In the calculation harmonic oscillator states for 
the occupied (hole) states and plane waves for the inter- 
mediate unbound particle states are assumed, incorpo- 
rating the correct energy and density dependence char- 
acteristic of a finite nucleus C/Fjv-matrix. In a similar 
way, one can construct the second-order self-energy con- 
tribution depicted in Fig. [TJc) which has an imaginary 
part below the Fermi energy and includes the coupling 
to onc-particle-two-hole (lp2h) states. 

Calculations of this kind require several basis transfor- 
mations, including the one from relative and center-of- 
mass momenta with corresponding orbital angular mo- 
menta to two-particle states with individual momenta 
and orbital angular momentum. Complete details can 
be found in Refs. [sHi- In practice the imaginary parts 
associated with diagrams (b) and (c) of Fig. [T] are em- 
ployed to obtain the corresponding real parts by employ- 
ing the appropriate dispersion relation. The resulting 
(irreducible) self-energy then reads 



^BHF + (F'-C ^2plh 



(6) 



i (Im S^pi;, 



in obvious notation. This self-energy is employed in the 
sp basis denoted by states \{k{£^)jmj]), characterized 
by wave vector, orbital, spin, total angular momentum 
and its projection (suppressing isospin). We note that 
the quantum numbers i,j and rrij are conserved and the 
self-energy does not depend on rrij . We will make use of 
isospin conservation for ^°Ca but include the Coulomb 
contribution when discussing results for protons. 



B. Solution of the Dyson equation 

The sp propagator in momentum space can be ob- 
tained from the following version of the Dyson equa- 
tion [ail 



G.,(fc,fc';i.;) = t_^G(°)(fc;i?) 



(7) 



+ G(°) (fc; £;)E^j (fc, k'; E)G'-°'^ (fc; E), 

where G''^\k;E) = {E — h'^k'^ /2m + irj)^^ corresponds to 
the free propagator and S^j is the reducible self-energy. 
The latter can be obtained by iterating the irreducible 
self-energy to all orders 

T.tj{k,k';E) ^j:}Jk,k';E) (8) 



h J dqq^Y,}^ {k, q- E)G^°^ [q; E)^,, (g, k'-E). 
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We first concentrate on energies below the Fermi energy 
that involve quantities that elucidate properties of the 
ground state. Below the Fermi energy the hole spectral 
function is determined by the imaginary part of the prop- 
agator 



1 



Sej{k;E) = -Ini Gij{k,k; E). 



(9) 



for the overlap functions is introduced. For an eigenstate 
of the Schrodinger-like equation [Eq. p^ ]. the so-called 
quasihole state labeled by aqh, the corresponding nor- 
malization or spectroscopic factor is given by [36j 



dE 



(17) 



For negative energies, the free propagator has no imagi- 
nary part and so according to Eq. ([7]) the spectral func- 
tion reads 

SiAk-E) = -G^°'>{k;E)lm Y.ej{k, k; E)G^"'> {k; E), 

TT 

(10) 

for energies where the imaginary part of the self-energy 
does not vanish. The total spectral strength at E for a 
given £j combination, 

poo 

Si,{E) = dkk^ Sej{k;E), (11) 
Jo 

yields the spectroscopic strength per unit of energy. 

The imaginary part of the CDBonn self-energy van- 
ishes between the maximum energy of lp2h and the min- 
imum energy of 2plh states. Inside this domain, discrete 
solutions to the Dyson equation are obtained from the 
standard version of the Dyson equation 



G,^ik,k';E)=^-^i^^G^'Hk;E) 



(12) 



+G^°Hk;E) / dqq^j:;^{k,q;E)Gi,iq,k';E). 



The solution for discrete poles utilizes the spectral rep- 
resentation of the propagator 



Gej{k,k';E) ^J:, 



E-{E. 



k'ij 



(13) 



where complete sets of states in the Azt 1 systems are 
inserted. The continuum solutions in the A ± I sys- 
tems are also implied in the completeness relations. The 
numerators of the particle and hole components of the 
propagator represent the products of overlap functions 
associated with adding or removing a nucleon from the 
A-body ground state. Following standard steps [s^ the 
eigenvalue problem corresponding to Eq. (jl2p reads 

where 



<^^^(fc)+ /dgg^S,*^.(fc,q;e-)0^.(<z)=£,;</.^_,.(fc), (14) 



E. 



A-l 



(15) 



for energies below the corresponding Fermi energy Sp 



E, 



The notation 



/5^</)^;,.(fc) = (vI/^Vfe^.l*o: 



(16) 



which is the discrete equivalent of Eq. pT|) . Discrete 
solutions in the domain where the self-energy has no 
imaginary part can therefore be obtained by expressing 
Eq. (|14p on a grid in momentum space and performing 
the corresponding matrix diagonalization. This also ap- 
plies to bound orbits above the corresponding Fermi en- 
ergy given by = E^^^ - Eq - We note that the DOM 
self-energy of Ref. [20[ contains an imaginary part below 
(hole domain) and above (particles) the average Fermi 
energy ep = {sp + ep)/2. 

The momentum distribution for a given £j is obtained 
from 



nij{k) = njj{k) + n'^^^^k) 



(18) 



where the continuum contribution is obtained by inte- 
grating the spectral function up to corresponding thresh- 
old 



dE Sej{k;E) 



(19) 



and the contribution of the discrete quasihole states reads 



(20) 



For protons the total momentum distribution (normal- 
ized by the number of protons Z) is obtained from 



(fc) = lE(2j + l)n,,(fc), 



(21) 



with a similar result for the neutron distribution. 

Information about natural orbits can be generated by 
determining the one-body density matrix. The contin- 
uum contribution to the one-body density matrix reads 



n^^(fc',fc) = i / dE Sej{k,k';E) (22) 



where 



Sej{k,k';E) ^ -G^°\k;E)lm Y.tj[k,k' ; E)G^°\k' ■ E) 

TT 

(23) 

corresponds to the nondiagonal spectral density. The 
one-body density matrix also receives a contribution from 
the quasiholes according to 



^{k',k) = Y^s^^4>7;{m,{n 



(24) 
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where S^j is the spectroscopic factor and (j}^j{k) are the 
quasihole eigen functions in momentum space. The total 
one-body density matrix is then given by 



nij{k', k) = nl^{k' , k) + n1j{k' , k). 



(25) 



By diagonahzing the one-body density matrix given in 
Eq. ((25|) one obtains the natural orbits for each ij combi- 
nation together with the corresponding occupation num- 
bers. It is therefore possible to write 



nij{k,k') = 



(26) 



with n"f-,(j)"f-{k) the corresponding occupation numbers 
and wave functions for natural orbit i. We note that 
these wave functions are normalized to unity. Results 
for quantities introduced in this section will be discussed 
in Sec. IIIII and wherever possible compared to the corre- 
sponding DOM quantities. 

In the language of many-body theory, the elastic 
nucleon-nucleus scattering is determined by the on-shell 
matrix element of the reducible self-energy YiiAk, k'; E), 
since it is directly related to the iS-matrix by [36| 



{ko\Sej{E)\ko)=e^'''^ 



(27) 



1 - 2TTi 



mka 



(fco|S«j(-B)|fco) 



where fco = \j2mE jh^ m is the nucleon mass, and E 
is the energy relative to the center-of-mass. The phase 
shift, 5ij, defined by Eq. (P7)) is in general a complex 
number. Its real part yields the usual phase shift and its 
imaginary part is associated with the inelasticity of the 
scattering process and denoted by 



_ „-2Im(5f,) 



(28) 



In general, the coupling to more complicated excitations 
in the self-energy implies a complex potential responsible 
for the loss of flux in the elastic channel, characterized 
by the inelasticities ijij . 

Because self-energy calculations at positive energy are 
rare, it is perhaps useful to include some relevant results 
in terms of the phase shifts 5ij for the quantities that will 
be discussed later . The scattering amplitude is given by 



4m7r^ , , 



ft2 



{k'm'Jj:{E)\km,), (29) 



with wave vectors of magnitude /cq. The matrix structure 
is usually represented by 



[f{0,^)]^T{e)i + cT-hg{e), 



(30) 



based on rotational invariance and parity conservation. 
The unit vector is given by n = A; x k' /\k x k'\, and cr is 
formed by the Pauli spin matrices. The relation between 
J' and g and the phase shifts determined by Eq. (I?fl) . 



can now be worked out yielding 

1}] P,(cos^) 



2ik 

n 



e=o 



(31) 



and 



sie) 



sm t 



E 



2iS, 



P'lico&e). (32) 



We employ the notation 5f^± = 5ij=i±^ and denotes 
the derivative of the Legcndre polynomial with respect 
to cosO. The unpolarized cross section reads 



da 



\Q\' 



(33) 



unpol 



Employing the partial wave expansions pip and p2l) and 
the orthogonality of the Legcndre polynomials, we find 



7rf. |(^+l){ 



l}+i{ 



1}| 



i=0 



-T 



2^+1 



+ l)|e' 



2iSr 



„2iS, 



2e + i 



(34) 



We can define partial elastic cross sections such that 



e=o 



which for a given £ read 



TT 
fc2 



+ 1) |e 



1 +e e 



2i5, 



1 



(35) 



(36) 



With complex potentials, and therefore complex phase 
shifts, it is possible to calculate the total reaction cross 
section 



E-- 



1=0 



with 



TT 



(37) 



(38) 



These results are derived by using the optical theorem 
that yields the total cross section from the imaginary 
part of the forward scattering amplitude [37| 



(39) 



The results presented here refer to neutron scattering 
when the self-energy has a finite range. 

Therefore, it is clear that at positive energies the prob- 
lem is completely reduced to solving the integral equation 
for the reducible self-energy given in Eq. ([S]). It should 
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be noted that the solution in momentum space automat- 
ically treats the non-locality of the reducible self-energy 
in coordinate space. In practice, the integral equation is 
solved in two steps. First the integral equation is solved 
by only including the principal value part of the nonin- 
teracting propagator. Subsequently, it is straightforward 
to employ the resulting reaction matrix to take into ac- 
count the contribution of the (5-function associated with 
the imaginary part of the noninter acting propagator. 

Previous calculations of the nuclcon self-energy for 
^^O 0-0 only included a limited number of sp partial 
waves. While this is sufficient for energies below the 
Fermi energy, results for differential cross sections at pos- 
itive energy quickly require many more partial waves with 
increasing energy. It is nevertheless possible to draw sen- 
sible conclusions by comparison with the corresponding 
partial wave contributions from DOM results that pro- 
vide an accurate fit to the experimental data. 



III. RESULTS 

A. Spectral Functions 

In Fig. [5] the ^3 /2 spectral function as a function of 
momentum is shown for different negative energies. The 
(^3/2 orbit is the last one that is mostly occupied in ^°Ca 
and the momentum content of the state near the Fermi 
energy exhibits no substantial strength beyond the naive 
expectation, similar to what was found for the Pi/2 orbit 
in ^^O before As ihustrated in Fig. the strength 
at higher momenta increases with decreasing energy, as 
expected and consistent with earlier calculations Q. We 
note that the presence of the high-momenta is not as 
pronounced as observed previously in ^^O in Refs. 0, 
In this earlier work the Bonn-B potential (ssj was em- 
ployed, whereas the present work employs the CDBonn 
interaction [10, [3l| which appears to be softer. We note 
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FIG. 2. Momentum-space spectral function for neutron 
quantum numbers at different energies in '"'Ca. 



10 



^ 10" 




!^ . J-" 



10" 



\ 10"^ 



10" 



10 



10" 





■4 •-' 



-150 



n ' 1 — 

100 -50 
E [MeV] 



FIG. 3. Comparison of spectral functions from present work 
(upper panel) with those from a DOM analysis using a non- 
local potential (lower panel). The results are for neutrons 
in ■'"Ca. The curves correspond to different ^j-combinations 
that are identified in the text. 



that the absence (presence) of high momenta near (far 
below) the Fermi energy is a simple consequence of si- 
multaneous energy and momentum conservation and is 
well-documented for nuclear matter [3l|. The behavior 
of the spectral functions in Fig. [5] is somewhat different 
from the DOM results discussed in Ref. 



22|. Whereas the 



microscopic strength exhibits a decrease of low momenta 
starting at -75 MeV, while the higher momenta exhibit 
similar strength, the corresponding DOM results show a 
more uniform behavior with little change in momentum 
profile as a function of energy. 

In the top panel of Fig. [3] the discrete and continuum 
contributions to the spectral strength [see Eq. (fTTj) ] are 
shown as a function of the energy for various ij channels 
for neutrons in ^°Ca. Results are shown for the follow- 
ing ij combinations, si/2 (solid), py,/2 (long-dash), pi/2 
(long-dash-dot), (dash), c?3/2 (dotted), f'j/2 (dash- 
dash-dot), and /5/2 (dash-dot-dot). In the bottom panel 
the corresponding DOM strength functions are shown for 
comparison. The DOM self-energy from Ref. [2^ has an 
imaginary part that ends at = ep and includes a sur- 
face term to account for long-range correlations (LRC). 
The resulting strength distribution is therefore continu- 
ous with sharp peaks near the Fermi energy. These peaks 
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are represented by discrete ones (normalized by spectro- 
scopic factors in the figure) for the CDBonn self-energy. 
Because the coupling to LRC is deemphasized, the imag- 
inary part ends at a much lower energy = —38 MeV. 
The DOM peaks closely correspond to the location of 
corresponding experimental sp states. The CDBonn self- 
energy generates energies for the peaks near the Fermi en- 
ergy that are not too different from the DOM result. For 
the deeply bound Si/2 peak, the CDBonn result underes- 
timates the binding by a substantial amount as compared 
to the DOM result which is consistent with a correspond- 
ing proton peak observed in the (e, e'p) reaction [40|. 

From Fig. [3l it is evident that the microscopic self- 
energy generates more sp strength at energies below the 
deeply-bound Si/2 peak than the DOM self-energy used 
in Ref. [2^ . Another important difference is that the sp 
strength decreases more quickly with increasing £ in the 
microscopic approach than in the DOM. We attribute 
this behavior to the non-locality of the imaginary part of 
the microscopic self-energy, whereas the imaginary part 
of the DOM self-energy is purely local. This feature has 
important consequences for the number of particles calcu- 
lated from the corresponding spectral functions according 
to 

N,aic = 5^(2.? + 1) f ' dE Sej{E). (40) 

Since the microscopic self-energy shows a substantial re- 
duction of sp strength with increasing the total particle 
number exhibits a good convergence with Ncaic = 19.3 
for neutrons when also the two partial waves with i ^ A 
are included. The corresponding DOM result exhibits 
a much slower convergence as illustrated in Fig. O and 
crosses 20 already at the orbit We therefore 

conclude that the non-locality of the imaginary part of 
the self-energy plays an important role in generating the 
correct number of particles from a DOM self-energy. 

B. Quasiholes and quasiparticles 

The quasihole energies are shown in Table H] They were 
obtained by diagonalizing the self-energy in momentum 
space while taking the energy dependence of the real part 
of the self-energy properly into account but disregarding 
the imaginary part. We also include the result for the 
static contribution of the self-energy labeled by BHF in 
the table. We note that the inclusion of the dispersive 
contribution moves the Osi/2 state up by almost 13 MeV. 
The cause of this huge shift must be attributed to the 
strong energy dependence of the diagram (c) in Fig. [T] 
which is very repulsive at the solution of the Osi/2 eigen- 
value. For other quasihole energies substantially smaller 
corrections of both signs are obtained. While not includ- 
ing the imaginary part of the self-energy does not yield 
the correct normalization for the Osi/2 state, we do find 
that the energy in the table is consistent with the loca- 
tion of the corresponding peak in Fig. |31 Results for the 



TABLE I. Quasihole energies for neutron orbits in "Ca. The 
second column shows the results from using the Hartree-Fock 
part only and third column the results from including the 
2plh and 2hlp terms in the self-energy. We also include DOM 
results and the position of the experimental sp levels near the 
Fermi energy. 

orbit BHF [MeV] Full [MeV] DOM [MeV] Exp. [MeV] 



OSi/2 


-56.1 


-43.6 


-56.1 




OP3/2 


-37.4 


-33.9 


-39.6 




Opi/2 


-34.7 


-31.7 


-34.9 




0d5/2 


-20.4 


-21.8 


-21.6 


-22.3 


lSl/2 


-18.1 


-19.6 


-17.4 


-18.3 


0d3/2 


-16.0 


-17.8 


-15.9 


-15.6 


0/7/2 


-4.3 


-7.1 


-9.8 


-8.4 


lP3/2 


-2.6 


-5.1 


-7.0 


-5.9 


lPl/2 


-1.2 


-3.5 


-5.4 


-4.2 



DOM self-energy employed in Ref. [22] are also included 
in the table as well as the experimental location of the sp 
orbits near the Fermi energy. The particle-hole gap of the 
CDBonn self-energy is more than 10 MeV, substantially 
larger than for the DOM at 6.1 MeV which is a little be- 
low the experimental result of 6.8 MeV. A common issue 
with microscopic self-energies is the underestimate of the 
spin-orbit splitting near the Fermi energy. For the split- 
ting of the d-states below the Fermi energy only 4 VIeV 
is obtained compared to 6.7 MeV found experimentally. 
The DOM generates 5.7 MeV for this quantity. Relativis- 
tic effects and core polarization [4]| or the importance of 
three-body forces |42l ] are usually invoked to repair this 
discrepancy but we note that only part of the sp strength 
resides in these levels, making the determination of the 
spin-orbit splitting more ambiguous. 

The spectroscopic factors identify the amount of sp 
strength residing near the Fermi energy and together 
with occupation numbers are shown in Table for the 
quasiparticle and quasihole states. The spectroscopic fac- 
tors are calculated according to Eq. ([TT]) after solving the 
Dyson equation (|14p . Occupation numbers for the quasi- 
particle or quasihole orbits arc obtained by double folding 
the momentum space wave functions with the one-body 
density matrix in momentum space given in Eq. (|25p . For 



TABLE II. Quasihole (quasiparticle) spectroscopic factors 
and occupation numbers for the CDBonn (CDB) self-energy 
compared to the corresponding DOM results. 



orbit 


ScDB 


Sdom 


ncDB 


riDOM 


lSl/2 


0.85 


0.66 


0.91 


0.88 


0^3/2 


0.87 


0.64 


0.92 


0.86 


0/7/2 


0.92 


0.67 


0.02 


0.11 


lP3/2 


0.93 


0.69 


0.02 


0.07 


lPl/2 


0.93 


0.73 


0.02 


0.06 
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the hole states near e p there is a reduction of sp strength 
of a httle more than 10% for the CDBonn calculation as- 
sociated mostly with the effect of SRC. For the particle 
states, the reduction of the sp strength corresponds to 
about 10%. It was observed in Ref. [1^ that the spectro- 
scopic factors calculated for orbits in the continuum are 
not reliable (and can be larger than 1) so these are not 
included in the table. All CDBonn spectroscopic factors 
are about 20% larger than in the DOM calculation. Since 
the DOM includes the coupling to low-lying excitations 
associated with collective effects of the nuclear surface, 
this simply reflects the important role of LRC that are 
responsible for this difference, and the resulting spectro- 
scopic factors are in good agreement with the analysis 
of the (e, e'p) reaction on '^''Ca ^^l- Differences in inte- 
grated strength as displayed by occupation numbers are 
less dramatic partially because the strength removed to 
lower energy from the quasihole (particle) peaks due to 
LRC in the DOM is recovered in the occupation num- 
bers [m. 
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FIG. 5. (Color online) CDBonn momentum distribution with 
imax ~ 3 (dashed) and Imax ~ 4 (solid) compared with the 
DOM result including all partial waves including f-j/2 (dash- 
dot) as obtained in Ref. [22l. The normalization of the curves 
is given by Atv J dk k^n{k) = 1 



C. Momentum distribution 

The total momentum distribution for neutrons in "^^Ca 
resulting from the CDBonn interaction and calculated ac- 
cording to the equivalent of Eq. (PT|) is shown in Fig.|3]by 
the solid line. All distributions arc multiplied by and 
normalized such that Air J dk k'^n{k) = 1. The quasihole 
contribution without the reduction due to spectroscopic 
factors is shown by the dashed line. Comparison of the 
correlated with this mean-field-like contribution (dashed) 
shows an appreciable presence of high-momentum com- 
ponents. A comparison with the DOM result (dash-dot) 
in Fig. [5l normalized as in Fig. |4l indicates that the 
DOM self-energy has slightly more high-momentum com- 



ponents than the microscopic one, although the latter 
has a larger quasihole contribution at high-momentum, 
as shown Fig.lB] The amount of strength above 1.4 fm~"i 
corresponds to 8% for the CDBonn calculation compared 
to 10% for the DOM. The convergence of the CDBonn 
result with orbital angular momentum is also illustrated 
in Fig. [5] exhibiting a satisfactory convergence when the 
(max = 3 result (dashed) is compared with the one for 
imax = 4 (solid). This convergence is less satisfactory 
for the DOM result [l^l which contains contributions 
from partial waves only up to and including the /7/2 
one. These results further confirm the importance of 
a non-local representation of the imaginary part of the 
self-energy which automatically leads to a more satisfac- 
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FIG. 4. (Color online) Momentum distribution for neutrons in 
*°Ca weighted by k^. Solid line represents the total momen- 
tum distribution including quasihole and continuum terms. 
The dashed line represents the quasihole result without re- 
ductions from spectroscopic factors. 
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FIG. 6. (Color online) Quasihole contribution to the mo- 
mentum distribution for the CDBonn (solid) compared to the 
DOM result (dashed). Both distributions are normalized ac- 
cording to 47r J dk k^n{k) = 1. 
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tory convergence with orbital angular momentum. As 
discussed in Sec. IIII A[ the energy dependence of the 
spectral function of the CDBonn potential in momentum 
space already suggests that it is a rather soft potential 
in comparison with the Bonn-B potential [38| that was 
employed for in Refs. 0-01 • While the spectroscopic 
factors for the aforementioned Bonn potentials in these 
nuclei arc similar, the CDBonn potential contains about 
4% of strength in the quasihole orbits above 1.4 fm~^, 
whereas for the Bonn-B potential this amount is much 
smaller. It would be interesting in future to investigate 
the corresponding behavior of a harder and local inter- 
action such as the Argonne V18 [43 |. 



D. Natural Orbits 
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FIG. 8. As in Fig. [7] but showing si/2 wave functions with 
one node. 



Results for natural orbits are obtained by diagonalizing 
the one-body density matrix which leads to Eq. pG]) con- 
taining the natural orbit functions in momentum space 
(shown in Figs. [TUHl) accompanied by the corresponding 
occupation numbers. In Tabic Hill results for the occupa- 
tion numbers of natural orbits are collected for the rele- 
vant partial waves. The results are quantitatively similar 
to previous results for ^^O @. More surprisingly, the 
results for the occupation numbers of the natural orbits 
that are predominantly occupied are very close to the 



DOM result reported in Ref. [22|. Only the orbits imme- 
diately adjacent to the Fermi energy exhibit more than 1 
or 2% deviation but never more than 8%. This is perhaps 
surprising in light of the rather large difference between 
the spectroscopic factors obtained with these different ap- 
proaches. Since both approaches are constrained by very 
different experimental data, it is therefore gratifying that 
the dominant occupation numbers for the natural orbits 




k[fm ] 



FIG. 7. (Color online) Comparison of wave functions for the 
51/2 quasihole result (solid) and the corresponding natural 
orbit (dashed) without a node. 



are so close. Differences of a few percent are observed for 
the smaller occupation numbers leading to a total occu- 
pation of these DOM orbits that is a few percent higher 
than in the case of the CDBonn. We ascribe this to the 
inclusion of LRC in the DOM. 

A comparison with natural orbits obtained for finite 
drops of ■^He atoms [1^ illustrates the substantial dif- 
ference between the underlying fermion-fermion interac- 
tions. The atom-atom interaction being much more re- 
pulsive and in a larger domain, relative to the volume per 
particle, leads to occupation numbers as small as 0.54 for 
the Is state in a drop of 70 '^He atoms, whereas the nu- 
clear interaction generates values close to 0.9 either in the 
DOM or from the CDBonn calculation. The difference is 
therefore mostly related to the much stronger repulsion 
between '^He atoms which e.g. in the liquid at satura- 
tion leads to a depletion of the Fermi sea of more than 
50% Nucleon-nucleon interactions typically gener- 

ate 10-15% depletion due to SRC [32|. 
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FIG. 9. As in Fig. [71 but with the corresponding ^3/2 wave 
functions. 
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TABLE III. Occupation numbers of natural orbits 



n 


Sl/2 


P3/2 


Pl/2 


^5/2 


d-i/2 


/7/2 


/5/2 


1 


0.882 


0.902 


0.898 


0.909 


0.919 


0.024 


0.025 


2 


0.910 


0.025 


0.025 


0.014 


0.015 


0.006 


0.007 


3 


0.016 


0.007 


0.008 


0.003 


0.004 


0.001 


0.001 


4 


0.004 


0.001 


0.002 


0.001 


0.001 


0.0003 


0.0005 


5 


0.001 


0.0004 


0.0004 


0.0001 


0.0002 


0.0001 


0.0001 


6 


0.0002 


0.0001 


0.0001 


0.0001 


0.0001 


< le-4 


< 1 e-4 




1.82 


0.94 


0.93 


0.93 


0.94 


0.03 


0.03 



Wave functions of natural orbits and quasihole states 
are compared in Figs.[7][Hl In Fig.[7]the momentum-space 
wave functions without a node for the Si^2 partial wave 
are shown for the deeply bound quasihole (solid) and the 
corresponding natural orbit (corresponding to an occu- 
pation number of 0.882). The natural orbit wave func- 
tion extends farther out in momentum space reflecting 
the contribution from lower energies that contain higher 
momenta. The quasihole wave function for the Si/2 state 
near the Fermi energy (solid) is compared to the natural 
orbit result (also with one node) in Fig. HI Again there is 
a substantial difference between the two wave functions 
although they both have the same occupation number. 
Contrary to the two si/2 results, the quasihole and 
natural orbit wave functions are essentially indistinguish- 
able, as shown in Fig. |9l A similar result was obtained 
(in coordinate space) for the DOM calculation of Ref. [22j 
and appears to be due to the presence of only one natu- 
ral orbit with a large occupation number unlike the Sjjj 
case. We note that in the DOM calculation of Ref. [23 
both quasihole si/2 wave functions were essentially iden- 
tical to the natural orbit results. 



E. Charge Distribution 

Although we focus mainly on neutron results for *°Ca, 
it is useful to study the charge density distribution ob- 
tained for the CDBonn potential. For this purpose, the 
Coulomb potential was incorporated into the calculations 
by first transforming the irreducible self-energy to coor- 
dinate space. The Coulomb potential from a uniformly 
charged sphere was then added, and a matrix inversion 
was performed to get the propagator in a similar way as 
in Ref. [22|. The radius of the sphere was taken to be 
Rc = 1.31^1/^, which is what was used in the DOM 
analysis. The point charge distribution was then ob- 
tained from 



(41) 



where Uij refers to the one-body density matrix in co- 
ordinate space. The resulting distribution was then cor- 
rected for the experimental charge distribution of a sin- 
gle proton and a single neutron as in Ref. (47| . The final 



charge distribution is shown by the dashed line in Fig.fTOl 
and compared to the experimental one obtained from the 
Fourier-Bessel analysis of Ref. [i^ . 

The mean square radius of the CDBonn self-energy is 
obtained from 



1 



(42) 



yielding a value of 3.29 fm compared to the experimen- 
tal result of 3.45 fm taken from Ref. [i^. We note that 
microscopic calculations usually underestimate the ex- 
perimental results (see e.g. Ref. [H for ^^O). Since LRC 
correlations are not adequately incorporated, it is pos- 
sible that an improved charge density is obtained when 
they are included. We note that a recent microscopic cal- 
culation of the matter density in ''"Ca also concentrates 
too much matter near the origin [i^ suggesting that it is 
useful to consider comparison with as many experimental 
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FIG. 10. Charge density distribution for ''"Ca from the CD- 
Bonn self-energy (dashed) compared to experiment (solid). 
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quantities as possible for a more detailed understanding 
of the quality of the many-body calculations. 

We note that the DOM charge density also contains 
too much charge near the origin [2^ even though LRC 
arc incorporated. The DOM self-energy was however 
constrained to reproduce the mean-square-radius of the 
charge distribution. It appears therefore essential to in- 
clude apart from SRC, as shown here for the CDBonn, 
both non-locality and LRC in the DOM to obtain a 
charge density that is in better agreement with the data. 
Indeed, future fits to experimental data in the DOM 
framework can utilize elastic electron scattering data di- 
rectly to fit the parameters for nuclei where such data 
are available. 



F. Energy of the ground state 

From the momentum distribution and the spectral 
function in fc-space, the neutron contribution to the 
ground state per neutron can be calculated by using: 

Ncalc 2Ncalc ^ J 2m 

+ + 1) / <^^^^ '^^ ^ Stj{k-E), 

calc J J —OO 

(43) 

where En is the total energy from the neutrons, and Ncaic 
[see Eq. (|40| ] is the calculated number of neutrons when 
all partial waves with i < A are included. With this limit 
on the number of partial waves Ncaic = 19.3 and a cor- 
responding energy per neutron of -8.25 MeV is obtained. 

As for the charge distribution calculation, the energy 
for the protons was generated in coordinate space. The 
density matrix and spectral function were then trans- 
formed back to the momentum space in order to use Eq. 
(j43|l with En replaced by Ep and Ncaic by the calculated 
proton number Zcaia which is found to be Zcaic = 19.5. 
The resulting energy per proton is -4.91 MeV, and the 
resulting total energy per particle is E{'^^Ca)/A = -6.56 
MeV. This result is 1.85 MeV per particle more attrac- 
tive than the DOM result in Ref. [22, but still more than 
2 MeV/ A higher than the experimental binding of -8.55 
MeV/ A. From the results of the spectral functions it 
is clear that more strength occurs in the continuum for 
the CDBonn at very negative energies than in the DOM 
calculation demonstrating the importance of these con- 
tinuum contributions to the total energy. 

This importance was also recognized in Ref. Q , where 
it was shown that about two-thirds of the binding is due 
to the continuum even though it represents only about 
10% of the particles. Since a different interaction is used 
and a heavier nucleus is considered in the present work, 
it is instructive to quantify the role of the continuum 
also for the CDBonn calculation. The quasihole contri- 
bution to the energy of the ground state are well sep- 



arated from the continuum except for the peak corre- 
sponding to the lowest Si/2 orbit as shown in Fig. [31 We 
therefore assessed the Sx/2 continuum contribution by in- 
tegrating the strength in Eq. ([321) up to -50 MeV. The 
total binding from the continuum contributions of all par- 
tial waves then amounts to -105.88 MeV compared to a 
total of -159.16, therefore representing 67% of the total, 
a very similar result to the ^^O calculations of Ref. [Bj. 
We note further that even though the CDBonn interac- 
tion is relatively soft, the cancellation between kinetic en- 
ergy (419.15 MeV) and potential energy (-578.32 MeV) 
is quite substantial. 

Recent calculations employing the unitary- mo del- 
operator approach generate -8.51/ A close to experiment 
for ^°Ca also using the CDBonn interaction [53|. The 
present calculation of -6.56 MeV /A is almost 2 MeV 
per particle less, pointing to the need of additional cor- 
relations and an improved treatment of the propaga- 
tors included in the present self-energy calculation. We 
note that only noninteracting propagators have been 
used in the construction of the self-energy and that self- 
consistency was not attempted. In addition, the proper 
treatment of LRC may also be relevant for certain ob- 
servables like the charge density. A microscopic treat- 
ment of LRC is available for the self-energy in which 
random-phase-approximation phonons are summed to all 
orders in a Faddeev technique and inserted into the self- 
energy [3l.[25l-|2i|. This FRPA method does not exphcitly 
include high-momentum components so a combination of 
the current method and the FRPA needs to be developed. 
Self-energy calculations also have the additional advan- 
tage that other observables of interest can be evaluated. 
Such observables are provided in particular by elastic- 
nucleon scattering data. We turn to an analysis of the 
corresponding results from the CDBonn self-energy in the 
following section. 



G. Neutron- Ca scattering 

In this section we make a direct comparison between 
CDBonn and DOM results, and experimental data for 
neutron- nucleus scattering, more specifically the n-|-*°Ca 
process. References to data can be found in Ref. [20l[2]|. 
The total cross section obtained from CDBonn self- 
energy is comparable to the experimental data for en- 
ergies between 5 to 10 MeV, while for larger energies 
it docs not reproduce the observations as shown by the 
solid line in Fig. 1111 However, currently our CDBonn 
self-energy only includes contributions up to imax = 4 
which is a reasonable limit for calculating properties be- 
low the Fermi energy. If we consider the DOM potential 
and only add the partial-wave contributions to the total 
cross section up to imax = 4, the corresponding result 
also fails to reproduce the experimental data, as shown 
by the dashed line in Fig. 1111 The most useful procedure 
is therefore to compare the DOM and CDBonn calcu- 
lation with the inclusion of the same number of partial 
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FIG. 11. (Color online) Comparison of the experimental cross 
section (data points) with the CDBonn calculation (solid line) 
for n+^^Ca. The dotted line corresponds to the DOM fitted 
cross section, and the dashed line represents the DOM result 
including up to ^ = 4 only. The cross section obtained from 
the CDBonn self-energy has contributions only up to ^ = 4. 



waves. It is then gratifying to observe that the micro- 
scopic calculation generates similar total cross sections 
as the DOM for energies where surface contributions are 
expected to be less relevant, i.e. above 70 MeV. Although 
calculations for more partial waves are in principle pos- 
sible, the amount of angular momentum recoupling cor- 
responding to the appropriate basis transformations be- 
comes increasingly cumbersome. Nevertheless, it is clear 



that the current limit of 



4 is insufficient to de- 



scribe the total cross section already at relatively low 
energy irrespective of whether surface effects (LRC) are 
properly included. 

To visualize in more detail the specific contribution 
from each partial wave, in Fig. [T^] we display cross sec- 
tions calculated up to a specific imax for CDBonn (upper 
panel) and DOM (lower panel). Despite the limitation of 
not having CDBonn values for £ > A, both potentials 
provide a comparable cross sections up to this angular 
momentum cut-off in the explored energy range. In the 
lower panel we also include the converged cross section 
of the DOM potential that provides a good description 
of the data. 

In Fig. [131 each panel displays the differential cross sec- 
tion as a function of the scattering angle for the indicated 
energies, covering 11, 26, 65, and 95 MeV for which data 
are available (shown by the dots) . The solid line in every 
panel corresponds to the CDBonn result and the others 
to the DOM, with the dashed line representing the con- 
verged cross section and the short-dashed line limiting 
the sum up to £max = 4. It is quite sobering to realize 
how far from the experimental data the microscopic cal- 
culations are, especially at the lower energies. This issue 
represents an important challenge to all many-body ap- 
proaches that are of interest to the study of rare isotopes, 
since only hadronic reactions are available to study these 
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FIG. 12. (Color online) Effect of adding subsequent partial 
waves for the total cross section from the CDBonn self-energy 
(upper panel) and corresponding DOM result (lower panel). 
The total DOM cross section with a converged contribution 
of the partial-wave decomposition is also shown as the solid 
line in the lower panel. 



and their analysis relies on presently unavailable optical 
potentials. At the higher energies (65 and 95 MeV) the 
differential cross section, while somewhat out of phase 
with the DOM result (up to £ = 4), is of similar magni- 
tude. Clearly, the missing details of LRC at lower energy 
are needed to get better agreement with the data, while 
the missing ^-contributions also become quickly notice- 
able. Future work clearly requires the calculation for 
higher ^-values, but should also focus on an improved in- 
clusion of surface dynamics in order to better reproduce 
the experimental nucleon-nucleus scattering data. 

In order to gain more insight into the similarities and 
differences between the CDBonn self-energy and DOM 
potentials, we also display the inelasticities [see Eq. ([28]) ]. 
The CDBonn's inelasticities rjgj for each partial wave are 
presented in Fig. [14] as a function of energy. As expected 
rjej is equal to 1 for energies below 25 MeV since the CD- 
Bonn potential is purely real from —38 to 25MeV. We ob- 
serve that the loss of flux above 25 MeV for each partial 
wave as a function of energy is larger for states than 
for £~ ones starting at energies ~ 50 MeV, which reflects 
the presence of a non-negligible spin-orbit content in the 
imaginary part of the CDBonn self-energy. As expected. 
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FIG. 13. (Color online) Differential cross section for indicated 
energies. Shown are the CDBonn result (solid) , DOM for the 
same imax = 4 (short-dashed), and the full DOM cross section 
(dashed). 
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FIG. 14. (Color online) In each panel, the inelasticity rjej 
for the indicated states as a function of the centre of mass 
energy. The solid (dashed) lines correspond to j — \£ + s\ 
~ \^ ~ s|) states. Below 25 MeV, the CDBonn potential 
is real, therefore the rjej are equal to one. Above 25 MeV 
the decrease of the r]ij indicates a loss of flux as the energy 
increases. 



this difference is more pronounced for higher orbital an- 
gular momenta. In the range of energies considered all 
the inelasticities are always larger than 0.7. 

We contrast these results with the inelasticities derived 
from the DOM potential for the same partial waves, as 
shown in Fig. 1151 The inelasticities derived from the 
DOM potential display little variation in the considered 
energy range, and saturate quickly and relatively close 
to 0.5. As anticipated the DOM potential is therefore 
more absorptive than the CDBonn self-energy. This dif- 
ference can be attributed to the lack of surface dynamics 
of the latter. We also note the absence of an appreciable 
imaginary spin-orbit potential in the DOM potential, in 
contrast with the CDBonn result. It is unclear whether 
the absence of an imaginary spin-orbit potential in the 
DOM is due to a compensating contribution that is not 
present in the current CDBonn self-energy. It is also un- 
clear without further analysis whether the apparent ab- 
sence of a need for a non-local potential in the analysis of 
elastic scattering is due to the truly local character of the 
self-energy or due to a conspiracy of other contributing 
factors. We perform some initial analysis into this issue 
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FIG. 15. (Color online) In each panel, the inelasticity rjij 
for the indicated states as a function of the centre of mass 
energy obtained from the DOM potential. The solid (dashed) 
lines correspond to j = \i + s\ {j = \l — s\) states. The 
short dotted line indicates 0.5, note that the DOM potential 
predicts inelasticities about 0.5 without significant variations 
with energy or orbital angular momentum. 



in the next section. 



H. Analysis of the CDBonn self-energy 

Improving the analysis of elastie scattering data above 
the Fermi energy and observables related to quantities 
below the Fermi energy in a DOM framework appears to 
depend sensitively on the treatment of non-locality in the 
imaginary part of the self-energy. It is therefore useful 
to gain some insight into the properties of microscopic 
self-energies which may offer guidance how to implement 
such features in the future. We therefore performed a few 
simple fits to represent the central part of the imaginary 
part of the CDBonn self-energy in coordinate space at a 
given energy assuming the following form of the potential 



WNL{r,r') = Wo^W)VW)H 



(44) 



We deviate from the standard Percy prescription for non- 
locality by employing square-root factors of the func- 
tion /(r) which is still represented by the conventional 



> 




r [fm] 



FIG. 16. (Color online) Diagonal part of CDBonn imagi- 
nary self-energy at 65 MeV (solid), and the corresponding 
parametrized self-energy (dashed). The results shown are for 
1^0. 



Woods-Saxon form factor 



fir) 



1 + e "0 



(45) 



with R ~ vqA^^^ . The function H determines the degree 
of non-locality and is assumed to be a gaussian following 
Ref. m 



H 



1 



^3/2^3 



exp 



(46) 



When the angular dependence in H is projected out, an 
analytic solution is obtained for each orbital angular mo- 
mentum i 



7ri/2/33 



X exp 



13' 



^'i~lfM^z), (47) 



where z — 2rr'//3^ and ji is a spherical Bessel function 
with a purely imaginary argument. The fact that an 
analytic projection is possible provided the motivation 
of the choice of Eq. (|33]). In arriving at the result of 
Eq. (|47| use has been made of the relation between the 
spherical Bessel functions and the Legendre polynomials 
Pf. 



3i{z) 



1 



dt e"^Pt{t). 



(48) 



We have chosen to fit the imaginary part at 65 MeV 
partly because we expect that only at such energies the 
imaginary part of the microscopic self-energy has real rel- 
evance since the role of LRC is expected to be diminished. 
In practise, this means that only the 1 = self-energy 
needs to be represented in terms of Eq. (|T7|) . If the choice 
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FIG. 17. (Color online) CDBonn imaginary self-energy at 65 
MeV integrated over r' (solid), and the corresponding inte- 
grated parametrized self-energy (dashed). The results shown 
are for £ — 0. 




FIG. 18. (Color online) Imaginary volume integrals for the 
CDBonn self-energy at 65 MeV (circles), and the correspond- 
ing result for the parametrized self-energy (squares). 



of Eq. is appropriate, the other ^-values will be ad- 
equately represented as well. In Fig. [16] we display the 
diagonal of the central imaginary part of the self-energy 
in coordinate space for € = for the CDBonn potential 
by the solid line. The fit according to Eq. pij) is quite 
satisfactory and given by the dashed line. Quantitative 
results for diffuseness aq, radius rp, and the non- locality 
parameter /3 are discussed below. 

Another useful check on the overall relevance of the 
parametrization of the non-local content of the poten- 
tial is to integrate over the variable r' in Eq. (|Tf|) to 
sample the nondiagonal components and compare with 
the corresponding integral for the CDBonn self-energy. 
The result of this procedure is identified by E*^^ and 
shown in Fig. [T7] as a function of r for the parametriza- 
tion (dashed) and CDBonn self-energy (solid) for orbital 
angular momentum £ ^ 0. This more stringent test in- 
cluding the sampling of non-diagonal components of the 
self-energy, still yields a satisfactory representation of the 
microscopic potential. It is interesting to note that the 
shape of the "local" potential is more reminiscent of a 
standard volume absorption. 

Another useful quantity to gauge the characteristic of 
an absorptive potential is the volume integral. For local 
potentials this quantity is well-constrained by experimen- 
tal cross sections [1^ [2l| . A recent analysis of the FRPA 
self-energy reveals useful insights into the comparison of 
microscopic self-energies with DOM potentials at lower 
energy [2J] . As in Ref. 2M '^'^ define the volume integral 
for a given orbital angular momentum i by 

J^{E) ^An J drr^ J dr'r'^ Im T.*i{r,r';E). (49) 

For a local potential it reduces to the standard definition 
of the volume integral. 

The implied ^-dependence of the chosen non-local po- 
tential leads to predictions for higher ^-values for this 



quantity once a fit to the £ ~ component of the self- 
energy has been made. The result of the corresponding 
volume integrals per nucleon are shown in Fig. [18] as a 
function of the £-values considered for the CDBonn self- 
energy. Wc employ dots for the CDBonn results and cir- 
cles for the predictions based on Eq. (|44]) . The agreement 
appears very satisfactory and may be useful to extract 
the properties of the CDBonn self-energy for even higher 
f-values without recourse to an explicit calculation. 

The properties of the imaginary part of the CDBonn 
self-energy in terms of its non-locality content are sum- 
marized in Table IIVI for four different energies, one be- 
low and three above the Fermi energy. In all cases a 
substantial imaginary part of the CDBonn self-energy is 
present at the chosen energies. The parameters are fit- 
ted at each energy to reproduce the essential properties 
of the self-energy including the volume integral for £ ~ 0, 
as discussed above for the case of 65 MeV. 

We observe that the values for the diffuseness generate 
standard values for the higher energies but are substan- 
tially larger for 49 and -76 MeV. The radius parameter 
is quite small below the Fermi energy but yields rather 
standard values at positive energy. The value of the non- 
locality parameter is quite a bit larger than typically as- 
sumed for real non-local potentials. Wave function cor- 
rections for non-locality in the analysis of (e, e'p) reac- 



TABLE IV. Parameters from non-local fits to the imaginary 
part of the proton self-energy at different energies. Wo is in 
MeV, ro,ao, /3 are in fm, and Jw is in units of MeV fm^. 



Energy 


Wo 


ro do 
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\Jw/A\ 


\Jw/A\ 


MeV 










CDBonn 


-76 


36.30 


0.90 0.90 


1.33 
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1.25 0.91 


1.43 


73 
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65 


13.21 


1.27 0.70 


1.29 
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FIG. 19. (Color online) Imaginary volume integrals for the 
CDBonn self-energy as a function of energy for different l- 
values: ^ = (solid), i = \ (dashed), i = 2 (short-dashed), 
^ = 3 (dash-dot), and ^ = 4 (dash- dash- dot). 

tions typically assume values of /3 = 0.85 fm [sij. The 
DOM analysis of Ref. [22| introduced a non-local Hartree- 
Fock potential to allow the calculation of additional prop- 
erties below the Fermi energy from the spectral functions 
that are the solutions of the Dyson equation. The ad- 
justed non-locality parameter in that work corresponded 
to 0.91 fm. 

We note that with increasing energy the non-locality 
parameter decreases suggesting a trend to a more local- 
ized potential. Since for a local potential there is no 
dcpendence of the volume integral, we have investigated 
the behavior of for different values in a wide energy 
domain. The results of this analysis are shown in Fig. 1191 
The degree of non-locality appears to be largest below the 
Fermi energy with a substantial separation between the 
different values. The result for ^ = also demonstrates 
that it is possible to have the "wrong" sign for the vol- 
ume integral. This can happen because the microscopic 
self-energy develops negative lobes off the diagonal and 
a positive volume integral cannot be guaranteed as a re- 
sult, as must be the case for a local potential. Although 
the imaginary part above the Fermi energy is negative, it 
is conventional to plot the in ragi nary volume integral as 
a positive function of energy [2l| . At positive energy the 
volume integrals for different i at first exhibit a spread 
although not as large as below the Fermi energy. Above 
300 MeV however, the curves apparently become similar 
suggesting a trend to a more local self-energy. 



IV. CONCLUSIONS 

The properties of the microscopic self-energy of nucle- 
ons derived from the realistic CDBonn interaction have 
been studied for ^°Ca. The calculation involves a two- 
step procedure starting with the calculation of a Gnm- 
matrix interaction in nuclear matter for a fixed energy 
and density. In a second step, the Fermi structure of the 



finite nucleus is incorporated by expanding the finite- 
nucleus G F AT-matrix in the nuclear matter one, including 
up to second-order terms. The self-energy is obtained 
by including the corresponding self-energy terms with 
imaginary parts above and below the Fermi energy, with 
associated real parts obtained from the appropriate dis- 
persion relations. The analysis of the solutions of the 
Dyson equation below the Fermi energy includes spectral 
functions calculated in momentum space, momentum dis- 
tributions, quasihole properties (including spectroscopic 
factors), natural-orbit properties, the nuclear charge den- 
sity, and the energy of the ground state of '"'Ca. 

An important motivation for the present work is to 
generate insight from microscopic calculations what func- 
tional forms of the nucleon self-energy can be employed 
fruitfully in the analysis of experimental data in the DOM 
framework. Recent DOM work has also focused on "'"Ca 
and a parallel comparison with DOM results is therefore 
pursued throughout the paper. 

Nucleon spectral functions for the CDBonn potential 
exhibit similar features as those from earlier work for 
^^O using the Bonn-B potential although the former in- 
teraction appears somewhat softer. This leads to a less 
pronounced presence of high-momentum components at 
very negative energies. The energy distribution of these 
momenta is somewhat different than the one generated 
by the DOM self-energy although the fraction of high- 
momentum particles is about 10% in both calculations. 

Since noninteracting intermediate states are employed 
in the CDBonn self-energy and therefore LRC are not 
well incorporated, there is no imaginary part in a sub- 
stantial region around the Fermi energy. As a result, only 
the lowest si/2 state is broadened in accordance with ex- 
periment, whereas all other quasihole states are repre- 
sented by discrete states. The DOM calculation exhibits 
a more realistic distribution of the sp strength including 
appropriate widths for p-states as well. The location of 
the quasihole states in the CDBonn calculation is in rea- 
sonable agreement with experiment but the particle-hole 
gap is larger than experiment. The associated spectro- 
scopic factors are close to 0.9 consistent with the 10% 
fraction of high- momentum nucleons. The DOM spectro- 
scopic factors are about 0.2 smaller since the DOM self- 
energy includes a strong coupling to the nuclear surface 
leading to better agreement with the analysis of (e, e'p) 
reactions. 

The calculation of natural orbits demonstrates that the 
largest occupation numbers are close to 0.9 very simi- 
lar to a recent DOM calculation even though substantial 
differences in spectroscopic factors occur, as discussed 
above. It appears that nuclear natural orbits always 
generate such occupation numbers in contrast with fi- 
nite drops of •^He atoms, where they can be substantially 
smaller in accordance with the much stronger repulsion 
of the underlying interaction. 

The nuclear charge density from the CDBonn self- 
energy exhibits too small a radius and too much charge 
at the origin but is otherwise not too dissimilar from the 
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DOM results. An important difference however, is the 
presence of a substantial non-local imaginary self-energy 
below the Fermi energy in the microscopic calculations. 
This leads to a good convergence with orbital angular 
momentum for the number of particles which amounts 
to 19.3 neutrons when imax = 4. No such convergence is 
obtained with DOM calculations on account of the local- 
ity of the imaginary self-energy, thereby overestimating 
the number of particles. We therefore conclude that the 
introduction of non-locality in the imaginary DOM po- 
tentials in the future is an essential ingredient. 

The distribution of high-momentum nuclcons from the 
CDBonn calculation leads to their large contribution of 
67% to the energy per particle in agreement with ear- 
lier observations for ^^O. The more realistic distribution 
of high-momenta leads to about 2 MeV more binding 
per nucleon than from the DOM self-energy but still 2 
MeV above the experimental result pointing to the need 
of an improved treatment of intermediate states in the 
self-energy and the consideration of higher-order contri- 
butions in the nuclear-matter fjjvj\/-niatrix interaction. 

We also performed an investigation of the CDBonn 
self-energy at positive energy. The lack of higher partial 
waves considered in the calculation leads to an under- 
estimate of the total cross section. A comparison with 
the same number of partial waves in the DOM calcula- 



tion demonstrates that above 70 MeV there is reason- 
able agreement between the two approaches suggesting 
that volume absorption is well represented by the present 
calculation but surface contributions (LRC) are missing. 
The differential cross section are not in good agreement 
with the corresponding DOM results (which generate 
good fits to the data when all contributing partial waves 
are included). This lack of agreement is worrisome but 
suggest that a lot more work is needed, in particular the 
inclusion of LRC. 

Finally, an analysis of the non-locality of the imag- 
inary part to the CDBonn self-energy reveals that its 
main properties can be quite well represented by a gaus- 
sian non-locality. Typical non-locality parameters are 
somewhat larger than those found in the literature. Vol- 
ume integrals indicate that non-locality is very impor- 
tant below the Fermi energy. Above the Fermi energy, it 
is initially substantial but appears to weaken at higher 
energies. 
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